UNIL-IPT-02-06 



Localizing gravity on a 't Hooft-Polyakov monopole in 

seven dimensions 



g ; Ewald Roesslf] and Mikhail Shaposhnikov 

O 

^ I Institute of Theoretical Physics 

^1 University of Lausanne 

CH-1015 Lausanne, Switzerland 

m 

> 
o 

cn 

\ We present regular solutions for a brane world scenario in the form of a 't Hooft- 

Polyakov monopole living in the three-dimensional spherical symmetric transverse space of 
a seven-dimensional spacetime. In contrast to the cases of a domain-wall in five dimensions 
and a string in six dimensions, there exist gravity-localizing solutions for both signs of 
QhI the bulk cosmological constant. A detailed discussion of the parameter space that leads 

^ ■ to localization of gravity is given. A point-like monopole limit is discussed. 

> 

X 



o 



Abstract 



^Email: ewald. roessl@ipt.unil.ch 

^ Email: mikhail . shaposhnikov@ipt . unil. ch 



1 Introduction 



Recently, there has been renewed interest in brane-models in which our world is represented as 
a 3 + 1- dimensional submanifold (a 3-brane) living in a higher-dimensional space-time |1|, 
This idea provides an alternative to Kaluza-Klein compactification and gives new insights 
to a construction of low energy effective theory of the fields of the standard model Q and 
gravity p, |^ . Moreover, it may shed light on gauge hierarchy problem p, and on cosmological 
constant problem 



-m- 



In string theory Standard Model fields are localized on D-branes - [|12|, whereas from the 



point of view of field theory brane model could be realized as a topological defect formed 
by scalar and gauge fields being a solution to the classical equations of motion of the coupled 
Einstein - Yang-Mills - scalar field equations. In the latter case one should be able to construct 
a solution leading to a regular geometry and localizing fields of different spins, including gravity. 

Quite a number of explicit solutions is already known. In five- dimensions a real scalar field 
forming a domain wall may serve as a model of 3-brane ||TB[. Higher dimension topological 
defects can be qualitatively different from a five-dimensional case from the point of view of 
localization of different fields on a brane. So, strings in six space-time dimensions were con- 



sidered in In particular, solutions corresponding to a thin local string together with 

fine-tuning relations (similar to the Randall- Sundrum domain-wall case) were found in (see 
also [20|). A numerical realization confirming the general results of |]18| in a singularity-free 
geometry and for the case of the Abelian Higgs model has been worked out in [0. Moving to 
even higher dimensions may be of interest because of a richer content of fermionic zero modes 
and because of more complicated structure of transverse space. In the framework of KK com- 
pactification monopoles in seven dimensions and instantons in eight dimensions were discussed 
in 



21, 22|, and brane-world scenarios in higher dimensions in [16], [23|-[27 



In [23] we considered a general point-like spherically symmetric topological defect as a 



model of 3-brane and formulated conditions that are necessary for gravity localization on it. 
A transition from a regular solution to the classical equations of motion to a point-like limit 
is in fact quite non-trivial for six and higher dimensions (see a detailed discussion for a string 
case in [p!9| , pO|). The aim of the present paper is to provide an existence proof of a possibility 
of gravity localization on a regular three-dimensional defect - 't Hooft-Polyakov monopole in 
seven dimensions, to study the parameter-space of a model that leads to gravity localization 
and to formulate exactly the meaning of point-like monopole limit. We confirm entirely the 
previous results, in particular, a possibility of gravity localization on a monopole embedded in 
a space with both signs of a bulk cosmological constant. 

The paper is organized as follows. In section |^ we present the 5*0(3) invariant Georgi- 
Glashow model (having 't Hooft-Polyakov monopoles as flat spacetime solutions ^8|,||2^) cou- 
pled to gravity in seven dimensions. The Einstein equations and the field equations are obtained 
in the case of a generalized 't Hooft-Polyakov ansatz for gauge and scalar fields. Boundary con- 
ditions are discussed in section |^, the asymptotic behavior of the solutions at the origin in 
section pTTI, at infinity in section p.2|. In section we give the relations between the brane 
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tension components necessary for warped compactification. Section ^ presents numerical results 
omitting all technical details. We give explicit sample solutions in section |4.1| and discuss their 
general dependence on the parameters of the model. We then present the fine-tuning surface 
(the relation between the independent parameters of the model necessary for gravity localiza- 
tion) in section L2. Sections |^ and ^ treat the Prasad-Sommerfield limit and the point-like 



monopole limit, respectively. While in the former case no gravity localizing solutions exist, in 
the latter case we demonstrate a possibility of the choice of the model parameters that leads to 
a fundamental Planck scale in TeV range and small modifications of the Newton's law, while 
well within the range of applicability of classical gravity. We conclude in section |^. In Ap- 
pendix A a derivation of the fine-tuning relations is given, whereas in Appendix B a discussion 
of the numerical details can be found. 



2 Field equations 

The action for the setup considered in this paper is a straightforward generalization of a gravi- 
tating 't Hooft-Polyakov monopole in 4 dimensions (which has been extensively studied in the 
past to the case of seven-dimensional spacetime: 

c _ c _l_ c fi ^ 

LJ gravity ^ ^ brane ' \ / 

Here Sg^^^ny is the seven-dimensional Einstein-Hilbert action: 



^gravity ^ J V i/^-"- ^5 J 



(2) 



g is the determinant of the metric gMN with signature ( — h + + + + +). We use the sign 



conventions for the Riemann tensor of Upper case latin indices M^N run over 0...6, 

lower case latin indices m, n over 4 ... 6 and greek indices /x, v over ... 3. The parameter M7 
denotes the fundamental gravity scale, A7 is the bulk cosmological constant and S^^^^^ is the 
action of Georgi and Glashow containing SU(2) gauge field W^j and a scalar triplet $° (we 
denote group indices by a, 6, c = 1 ... 3) : 

S,r.n. = J d\^gCm with £^ = -Ic'i.^G^*^^ - - ^ {^^^'^ - rff , 

(3) 

where rj is the vacuum expectation value of the scalar field and Vm is a covariant derivative, 

Vm^~^ = 9m<J>" + e e""^'' W~^<^'. (4) 

Furthermore one has 

Gl,j, = OmW'^ - dr,W!, + e e^' w'^W'^. (5) 
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The S0{3) symmetry is spontaneously broken down to U{1). The monopole corresponds to 
the simplest topologically nontrivial field configuration with unit winding number. The Higgs 
mass is given by ttih = rj V2X. Two of the gauge fields acquire a mass mw = erj. 



The general coupled system of Einsteins equations and the equations of motion for the scalar 
field and the gauge field following from the above action are 

Rmn — -QmnR + -j^Qmn = j^Tmn , (6) 
-^Vm{V^V^^<^'^) = A - 77^) , (7) 

^^^^Af(v^G™) = -ee''''(v''^'^^\ (8) 

where the stress-energy tensor Tmn is given by 

2 S S 

Tmn = ~ ^ — - ^gMN ~ ^ML^N + T'm^'^T^N^"^ + dMN^m ■ (9) 

We are interested in static monopole- like solutions to the set of equations respecting 
both, 4i5-Poincare invariance on the brane and rotational invariance in the transverse space. 
The fields and W^.^ (and as a result Gmat) should not depend on coordinates on the brane 
x'^. The brane is supposed to be located at the center of the magnetic monopole. A general 
non-factorisable ansatz for the metric satisfying the above conditions is 

ds^ = M\p)g^^Jdx''dx'' + dp^ + L\p) {de^ + sin^ 6 d^^) , (10) 

where gj^J is the four-dimensional metric that satisfies the 4D Einstein equations with an 
arbitrary cosmological constant Apj^yg In this paper we will only consider the case of 

Apj^yg = and we take gj^J to be the Minkowski metric rj^^ with signature ( — h -I--I-). 
With the use of spherical coordinates for transverse space 

= [sin 6 cos if, sin 6 sirup, cos 6) , (11) 
= (cos^^ cos cos 6* sin 9?, — sin ^) , (12) 
e'^ = (— siny9, COSV9, 0) , (13) 

the 't Hooft-Polyakov ansatz is : 

W, = 0, We = -'-^^e-, H>, = i^^sin^e^ W, = $=^e-, (14) 

e e ep 

with fiashes indicating vectors in internal SO (3) space. 
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Using the ansatz (|I^) for the fields together with the metric (|TD|) in the coupled system of 
differential equations (|D-(|D gives 



M" C" 



M' a 




M'^ 


+ 
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MC 


+ 3 


M2 






M' C 
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M'^ 
M2 
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1? 
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~ £2 



/2 
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4 \ h4 h6 
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J" + 2 
K" + 



M' 
H 



£2 



/:r £2 
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(15) 
(16) 
(17) 

(18) 
(19) 



where primes denote derivatives with respect to the transverse radial coordinate r, rescaled by 
the mass of the gauge boson mw'- 

r = mwP = vep. (20) 
All quantities appearing in the above equations are dimensionless, including a, (3 and 7: 



with 



" " ' mw) ' ^ Ml ' ^ e2?74 ' 627^4 



n = 6^Jo. TP = f,, T^ = f9, T^ = U, C = Lmw, J 



H 

Tje p 



The dimensionless diagonal elements of the stress-energy tensor are given by 



£2 



2CF 



+^J'V^+^(J^-1)^ 



£2 



K 



/2 



:i-ir2)2 



+ lj'2_:^_"(^2_i)^^ 



£2 



2£4 ' 2" £2 

« / ,0 . N 2 
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'^-^'^'-ij'^-^(j^-l)\ 



2£4 



(21) 
(22) 

(23) 

(24) 
(25) 



The rotational symmetry in transverse space implies that the (6*^) and the (v^V^) components 
of the Einstein equations are identical (and that = eo). Equations ([TSD - (|17D are not 
functionally independent [^. They are related by the Bianchi identities (or equivalently by 
conservation of stress-energy Vm^jv^ = 0). Following the lines of p3| we can define various 
components of the brane tension per unit length by 



Pi = - I drMirYCirfeiir) . 




(26) 
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The Ricci scalar and the curvature invariants i?, RabR^^, Rabcd R^^"^^ ^ CabcdC^^'-^^ 
with Cabcd being the Weyl tensor are given by 



R 2 2C'^ 16 CM' 12 M'^ AC" 8 M" 



ml, C? C? CM M2 C M ' 

RabR^^ _ 2 AC'^ 2C''^ IQC'M' 16 C'^M' ASC'^M'^ 

mtv ~ C'lC^^lC^ CM ^ C^M^ ^ 

ASC'M'^ 36 M'^ AC" AC'^C" 16 C'M'C" 6 C"^ 

+ —TT-. ^ + — + ^TT^T + -^z:r- + 



(27) 



16 C'M'M" 2AM'^M" 16 C" M" 20 M"^ 

+ + £M ^ M2 ' ^^^^ 

RabcdR^^^^ _ A _8C'^ AC'^ 32 C'^M'^ 2AM'^ 8 C"^ 16 M"^ 

Cabcd C^^^^ 8 16 8 ^ 128 £' M' 128 C'^M' 16 M'^ 



3£4 3£4 3£4 i5/:37Vf IS^^M 5£2M2 

208 r^M'^ 6AC'M'^ 2AM'^ 32 C" 32 C'^ C" 
15 £2 M2 5/:M3 ^ 5M4 ^ 15 £3 15£3 

64/:'M'/:" 32 M'^C" 6AC"^ 32 M" 32 C'^ M" 
15 £2 M ^ 5CM^ ^ 15 £2 " 15 £2 m + 15 £2 m + 

64 /:'M'M" 32 M'^M" 128 /:"M" 64 M"^ 
15 C M2 5M3 15 /:M ^ 15M2 ■ ^ ' 

They must be finite continuous functions for regular geometries we are interested in. 



3 Boundary conditions and asymptotics of the solutions 

The boundary conditions should lead to a regular solution at the origin. Thus we have to 
impose 

M\r=o = l, M'\r=0 = 0, C\r=0 = 0, £'|.=0 = 1, (31) 

for the components of the metric, where the value +1 for M\r=o is a convenient choice that can 
be obtained by rescaling of the brane coordinates. 

What concerns the gauge and the scalar fields, the boundary conditions for them are the 
same as for a monopole solution in the fiat space-time, ||28|| : 



J(0) = 0, limJ(r) = l, (32) 
K{0) = 1 , lim K{r) = 0. (33) 
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Finally, a requirement of gravity localization reads 

oo 

j M{rfC{rfdr = M|, < oo , (34) 



oo 



what puts a constraint on the behavior of the metric at infinity. 

The magnetic charge of the field configuration can either be determined by comparing the 
stress-energy tensor with the general expression given in or by a direct calculation of the 
magnetic field strength tensor p8[|: 



, $ ■ Gmn 1 

/MAf — 



$ ■ {Vm^ X Vn^^ . (35) 



The only nonzero component of Qmn is Qoip = Either way gives Q = ^. 



3.1 Behavior at the origin 

Once boundary conditions at the center of the defect are imposed for the fields and the metric, 
the system of equations (|TB|)-(|Ty|) can be solved in the vicinity of the origin by developing the 
fields and the metric into a power series in the (reduced) transverse radial variable r. For the 
given system this can be done up to any desired order. We give the power series up to third 



order in r 






M(r) 


= 1 - Ir^ /? (a + 4 7 - 6 K"iof) + 0{r') , 


(36) 


£(r) 


= r + —r^(3 (a + 4^-30J'{0f-66K"{0f) + 0{r^) , 
360 V / 


(37) 


J(r) 


-rJ'(O) 1 .3j,(o)-9« + «/3 + 4/37 + 6/3J'(0)Vl8ir"(0) + 6/5ir'W| ^.p^ 

90 


,(38) 


K{r) 


= l + ir2ir"(0) + O(r^). 


(39) 



It can easily be shown that the power series of M(r) and K{r) only involve even powers 
of r whereas those of C{r) and J(r) involve only odd ones. The expressions for C{r) and J(r) 
are therefore valid up to 5*^^ order. One observes that the solutions satisfying the boundary 
conditions at the origin can be parametrized by five parameters {a, J'{0), K"{0)). For 
arbitrary combinations of these parameters the corresponding metric solution will not satisfy 
the boundary conditions at infinity. Therefore the task is to find those parameter combinations 
for which ( ^4]) is finite. For completeness, we give the zero-th order of the power series solutions 
for the stress-energy tensor components and the curvature invariants at the origin: 

eo|r=o = -^(a + 6J'(0)2 + 6K"(0)2) , (40) 

e,\r=o = ee\r=o = {a + 2 J'{Of - 2 K"{Of) . (41) 
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R (3 



^ (7 a + 28 7 + 30/(0)' + 18 /^"(O)') , (42) 



Rab R P 



100 



(^7 + 56 a 7 + 112 7^ + 60 a J'{of + 240 7 J'{of + 300 /(O)^ 
+ 12 (^3 a + 12 7 + 70 J'(0)') K"{of + 732 ir"(0)^) , (43) 



J^ABCD P 



mf^ 300 



(^17 + 136 a 7 + 272 7^ - 60 a J'(0)' - 240 7 J'(0)' + 900 J'(0)^ 
-36 (^9 a + 36 7 - 110 J'(0)') fs:"(0)' + 4932 K"{<df^ , (44) 
'-^SE2£!^ - 40^-6 (.'(Of. 3A-W))\ (45, 

3.2 Behavior at infinity 

The asymptotics of the metric functions M(r) and C{r) far away from the monopole are PB| : 

M(r) = Mo and /:(r) = Cq = const , (46) 

where only positive values of c lead to gravity localization. This induces the following asymp- 
totics for the stress-energy components and the various curvature invariants: 

lim eo(r) = lim ep(r) = --^ , (47) 



r— >oo 



lim eg (r) = . (4^ 



lim^^ = -5c2 + 4' (49) 



hm J = 5c + -4 , (50) 

r^oo LI 

RABCD{r)R^^^^ = 5^ 4 

i^W ^ ''-0 



hm """"^ ^ = -0^ + -;^, (51) 



j.^ a4BCD(r)C;^f^ ^ ^ + (52) 

771^ 6 3 /^o /^o 

The parameters c and £0 are determined by Einsteins equations for large r and are given by 



1 5 I . I 32 



Only the positive signs of the roots lead to solutions with both > and 7^ > 



In order to obtain some information about the asymptotic behavior of the fields J(r) and 
K{r) at infinity we insert relations (|46|) into (18) and (19). Furthermore we use K{r) = 6K{r) 
and J(r) = 1 — 6J{r) with 6K -C 1 and 5 J ^ 1 for 1 <^ r to linearize these equations: 



6J" -2c6J' -2a6J 
6K"-2c6K'+ 

By using the ansatz 6K = Ae~^^' and 6 J = Be"^"^ we find 



0, 
0. 



(55) 
(56) 



/Ci 2 = -C± - 



£2 



jl,2 = -C ± Vc^ + 2Q! 



(57) 



To satisfy the boundary conditions we obviously have to impose /c > and j > 0. We distinguish 
two cases: 



1. c > Gravity localizing solutions. In this case there is a unique k for < 1 with the 



positive sign in (p7|). 



2. I c < I Solutions that do not localize gravity. 

Both solutions of (^) are positive. 
In this case there are no real solutions. 
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There is a unique solution with the positive sign in 
The important fc-value here is k = —2c. 



It can be easily shown that for large enough r the linear approximation (|56| ) to the equation of 
motion of the gauge field is always valid. Linearizing the equation for the scalar field however 
breaks down for 2k < j due to the presence of the term oc K^J in ([T8|) . In that case J(r) 
approaches 1 as e'"^^^. For solutions with c > (the case of predominant interest) a detailed 
discussion of the validity of 2k > j gives: 



1. I g = I Prasad-Sommerfield limit One has 2k > j. The asymptotics of the scalar 
field is governed by e~^'^. 



2. < a < 2 The validity of 2k > j depends on different inequalities between a, -k and 



a + -k<l 



2k >j . 
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(b) 



a+ -k>l 

'-0 



and 



In this case 2k > j is equivalent to < 



-I a i 



l-a- 



where equality in one of the equations implies equality in the other. 

2k < j and the scalar field asymptotics is governed by e"^'^''. 



3. a>2 



2k < j. Also in this case the gauge field K{r) determines the asymptotics of 



the scalar field J{r). 



3.3 Fine-tuning relations 

It is possible to derive analytic relations between the different components of the brane tensions 
valid for gravity localizing solutions. Integrating linear combinations of Einsteins equations 
(|T5|)-(|l7|) between and oo after multiplication with M'^(r)£^(r) gives: 



1^0 - N 



(3 



oo oo 

j M^dr = y" (1 - K') 



dr , 



(5^ 



27 / M^C'^dr 





00 



a 



j (1- J2) M^Cdr. 



(59) 
(60) 



To obtain the above relations integration by parts was used where the boundary terms dropped 
due to the boundary conditions given above. A detailed derivation of these relations is given 
in the appendix A. 
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4 Numerical solutions 



Details of numerical integrations are given in the Appendix P; in this section we will discuss 
the results only. 

4.1 Examples of numerical solutions 

Fig. |I| and Fig. ^ show two different numerical solutions corresponding to positive and neg- 
ative bulk cosmological constant, respectively. Both solutions localize gravity c > 0. The 
corresponding parameter values are (a = 1.0000000, (3 = 5.50000000, 7 = —0.05434431) and 
(a = 1.0000000, p = 3.50000000, 7 = 0.02678351). Figs. | and | and Figs. | and | show 
the corresponding components of the stress-energy tensor and the curvature invariants. For 
high values of (3 the metric function C{r) develops a maximum before attaining its bound- 
ary value. Gravity dominates and the volume of the transverse space stays finite. For lower 
values of f3 the transverse space has infinite volume and gravity can not be localized. See 
Figs. 1^ and || for a solution that does not localize gravity (c < 0) and corresponds to 
(a = 1.0000000, p = 1.80000000, 7 = 0.04053600). 
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Figure 1: Gravity- localizing solution with negative bulk cosmological constant corresponding to 
the parameter values {a = 1.00000000, (3 = 5.50000000, 7 = -0.05434431) 
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Figure 2: Stress-energy components for the solution given in Fig. [|. 
(a = 1.00000000, 13 = 5.50000000, 7 = -0.05434431) 
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Figure 3: Curvature invariants for the solution given in Fig. 
(a = 1.00000000, /3 = 5.50000000, 7 = -0.05434431) 
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Figure 6: Curvature invariants for the solution given in Fig. 0. 
(a = 1.00000000, p = 3.50000000, 7 = 0.02678351) 
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Figure 7: Example of a solution that does not localize gravity (c < 0) corresponding to the 
parameter values {a = 1.00000000, (3 = 1.80000000, 7 = 0.04053600) 
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Figure 8: Stress-energy components for the solution given in Fig. 0. 
(a = 1.00000000, p = 1.80000000, 7 = 0.04053600) 

4.2 The fine-tuning surface 

Fig. H shows the fine-tuning surface in parameter-space. A point on this surface (oq, Po, 70) 
corresponds to a particular solution with the metric asymptotics (^) for both values of the 
sign of c. The bold line separates gravity localizing solutions c > from solutions that do 
not localize gravity c < 0. The parameter space has been thoroughly exploited within the 
rectangles Aa x Ap = [0, 10] x [0.5, 4.0] and Aa x A/3 = [0, 1] x [3.5, 10] in the (a, /3)-plane. 
The series of solutions presented in section |4.1j can be used to illustrate their dependence on 



the parameter (3 (strength of gravity) for a for a fixed value of a = 1. 

It can be seen from Fig. ^ that for every fixed a there is a particular value of (3 such that 
c equals zero, which is the case for all points on the solid line shown in Fig. ^. By looking at 
eq. (0) we immediately see that c = is equivalent to /9^7 = -^t^tu = |- We will discuss the 
limit c -C 1 in more detail in section ^ It will turn out to be the most physical case where 
the monopole can be considered to be point-like since the fields attain their vacuum values 
much earlier than the metric goes to zero outside the core. A solution corresponding to that 



case is given in Fig. |T^, where gravity (parametrized by (3) is just strong enough to provide a 
finite volume for transverse space. If for fixed value of a we increase (3 (starting from c = 0), 
c becomes more and more positive, the Planck mass becomes smaller and the monopole size 
increases, see Figs. |1] and ^. If on the other hand /? is decreased (from c = on), c becomes 
more and more negative and the metric M(r) blows up exponentially, as in Fig. |^. Gravity is 
no longer strong enough to provide for a finite Planck mass. 

For a=0 (the Prasad-Sommerfield limit, to be discussed in section H) it can be read off from 
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Figure 9: Fine-tuning surface for solutions with the metric asymptotics ^^). The bold line 
separates solutions that localize gravity (c > 0) from those that do not (c < 0). Numerically 
obtained values 0/7 = are plotted as a function of a = and (3 = f^. 

Fig. ^ that there are no solutions that locahze gravity. All a = solutions lie in the c < part 
of the surface. 



5 Prasad-Sommerfield limit (a = 0) 

The Prasad-Sommerfield limit {a = 0) was exploited numerically for /5-values ranging from 0.4 
to about 70. The corresponding intersection of the fine-tuning surface Fig. ^ and the plane 
a = is given in Fig. There exist no gravity localizing solutions as the separating line in 
Fig. ^ indicates. The point in Fig. |1^ corresponds to the solution shown in Figs. ^ and [12. 
One sees that 7 tends to zero for (3 going to infinity and that 7 tends to —00 for (3 going to 
zero. 
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Figure 10: Section of the fine-tuning surface Fig. corresponding to the Prasad-Sommerfield 
limit a = 0. None of the shown combinations of (3 and j-values correspond to solutions that 
lead to warped compactification. The point indicates the sample solution given below. 
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Figure 11: Sample solution for the Prasad-Sommerfield limit corresponding to the parameter 
values {a = 0.00000000, /3 = 3.60000000, 7 = 0.02040333) 
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Figure 12: Stress-energy components for the solution given in Fig. |I7|. 
(a = 0.00000000, (3 = 3.60000000, 7 = 0.02040333) 

6 A point-like monopole limit - physical requirements 
on solutions 

The fine monopole can be characterized by c ^ 1. As anticipated in section |4.2| , c = in 
eq. (|53|) immediately leads to /5^7 = ^2^w = |- Hence we deduce that the fine- monopole limit 
can not be realized for a negative bulk cosmological constant. In addition, from it follows 
that Co = The c <^ 1 limit is qualitatively different from its analogue in the 6D-string case 
[ p^ . The solutions do not correspond to strictly local defects. The Einstein equations never 
decouple from the field equations. Due to the particular metric asymptotics (^Bp, the stress- 
energy tensor components tend to constants at infinity in transverse space. In the GD-string 
case the fine-string limit was realized as a strictly local defect having stress energy vanishing 
exponentially outside the string core. Despite this difference the discussions of the physical 
requirements are very similar. In the following we show that in the fine monopole-limit the 
dimension-full parameters of the system A7, Mj, mw, A, e can be chosen in such a way that all 
of the following physical requirements are simultaneously satisfied: 

1. Mj, equals (1.22 • lO^^GeV)^. 

2. The corrections to Newtons law do not contradict latest measurements. 

3. Classical gravity is applicable in the bulk. 

4. Classical gravity is applicable in the monopole core (r = 0). 
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To find solutions with the above mentioned properties it is possible to restrict oneself to a 
particular value of a, e.g. a = \. This choice corresponds to equal vector and Higgs masses 
my/ = fTtH- Even though extra dimensions are infinite, the fact that M(r) decreases exponen- 
tially permits the definition of an effective "size" tq of the extra dimensions: 

M = Moe"5'' = Moe ^ ^ _ (61) 

cmw 

In order to solve the hierarchy problem in similar lines to P] we parametrize the fundamental 
gravity scale as follows: 

M7 = /€lO^GeV. (62) 
n = 1 then sets the fundamental scale equal to the electroweak scale. 



1. The expression for the square of the Planck mass M|, can be approximated in the fine- 
monopole limit by using the asymptotics (|^) for the metric in the integral (|3^) rather 
than the exact (numerical) solutions. This gives 

M|«1^M.^£^1. (63) 

By using one of Einstein's equations at infinity 

^0 = 77^^ ' (64) 



4c2 + 2/^7 

and by developing to lowest order in c one finds 

Numerical solutions for c — and a = \ converge to the following approximate parameter 
values 

(3 = 3.266281 , 
7 = 0.0468665 , 
= 0.4999995 , 
Mo = 0.959721 . (66) 

The above values were obtained by extrapolation of solutions to c = 0, see Figs. |TB| to |T^. 
Therefore the relative errors are of order 10^^ which is considerably higher than average, 
relative errors from the integration which were at least 10~^. 

Note that jS'^j tends to | in the fine-monopole limit (see Fig. |l^). Neglecting all orders 
different from ^ in the above expression for Mp and using (pl| ) one has 

5/2 

= 1.1- 10-5 ^GeV, (67) 
where Tq has been parametrized by ro = "'^g""" ~ 10^^ GeV^^ ■ |. 
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Figure 13: Behavior of (3 in the fine-monopole limit c — > for a 



7=^ 



0.02 0.04 0.06 O.C 



Figure 14: Behavior of 7 in the fine-monopole limit c ^ for a 



2 ■ 



Figure 15: Behavior of jS'^'-f 



^T^w in the fine-monopole limit c —> for a 
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Figure 16: Behavior of Mq in the fine-monopole limit c ^ for a = \. 



2. Since Newtons law is established down to 0.2 mm [^] we simply need to have ^ > 1. 

3. In order to have classical gravity applicable in the bulk we require that curvature at 
infinity is negligible with respect to the corresponding power of the fundamental scale. 
By looking at eqs. (|49|)-(|52D we see that we have to impose 



(? <^ My and 



w 



r2 



(68) 



Using ( |64D and the first of the above relations, the second one can immediately be trans- 
formed into 



One then finds 



p / 3 \ 1/2 

2 ■ 10"^^- < 1 and 4 ■ lO'^ ( — 



«1, 



(69) 



(70) 



where /3 and 7 have again been replaced by their limiting values (|66|) for c — 0. 

4. Classical gravity is applicable in the monopole core whenever the curvature invariants i?^, 
RabR^^, RabcdR^^^^ and CabcdCabcd are small compared to the forth power of the 
fundamental gravity scale. Since these quantities are of the order of the mass (see 
we have 



Using (|62D and (q^) one finds 



1.1.10-8 I 1- 



3\ 1/2 



< 1. 



(71) 



(72) 
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Finally, the fine-monopole limit requires 



< 1. 



ro mw 



(73) 



Again (pTQ implies 



't3\ 1/2 

1.8-10-M^l <i 



(74) 



It is now easy to see that for a wide range of parameter combinations all these requirements on 
^ and K, can simultaneously be satisfied. One possible choice is ^ = 100 and k = 1. This shows 
that already in the case « = | there are physical solutions corresponding to a fine-monopole 
in the sense of eq. (^) respecting all of the above requirements ([!])-(§)• Fig. |l^ shows the 
fine-monopole solution corresponding to the lowest c-value in the sequence shown in Figs. |T^ 
to M. 



1.75 - 



L (r) 



1.5- 



1 .25 



J(r) 



. 75 



0.5 - 



0.25 - 




M (r) 



10 



15 



20 



Figure 17: Gravity-localizing solution in the fine-monopole limit for c = 1.728 ■ 10 ^ and corre- 
sponding parameter values of {a = 0.50000000, (3 = 3.27000000, 7 = 0.04676000) 



7 Conclusion 

We have demonstrated in this paper that it is possible to generalize the idea of warped com- 
pactification on a topological defect in a higher dimensional spacetime to n = 3 transverse 
dimensions by considering a specific field theoretical model. Numerical solutions were found 
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for the case of a monopole realized as a 't Hooft-Polyakov monopole. This generahzation turned 
out to be non-trivial since (at least) in the rotational invariant transverse space setup consid- 
ered, Einstein's equations don't seem to admit strictly local defect solutions. Even though the 
transverse space still approaches a constant curvature space at spatial infinity, it is not neces- 
sarily an anti-de-Sitter space. Both signs of the bulk cosmological constant are possible in order 
to localize gravity. We considered a fine monopole limit in the case a = 1/2 {rriyy = ttih) and 
verified that the model proposed is not in conflict with Newtons law, that it leads to a possible 
solution of the hierarchy problem and that classical gravity is applicable in the bulk and in the 
core of the defect. Even though stability should be guaranteed by topology, small perturbations 
around the monopole background should be considered as well as quadratic corrections to the 
Einstein-Hilbert action. 

Acknowledgments: We thank T. Gherghetta, H. Meyer, S. Randjbar-Daemi, P. Tinyakov, 
S. Wolf and K. Zuleta for helpful discussions. This work was supported by the FNRS grant 
20-64859.01. 



A Derivation of Fine-tuning relations 



By taking linear combinations of Einstein equations (15) to (|T^ one can easily derive the 
following relations: 

M{rYC{ry = -5(27 + 60-6,-26,), (75) 
M{rYC\r) ZM^ = --(27-4eo-e, + 3e,). (76) 

By multiplying with M (r)^£(r)^, integrating from to 00 and using the definition of the brane 
tensions (E6D we obtain 



M{rfM\r)C{r)X = T M(r)^/:(r)2rfr + ^ (/io - /ip - 2/i,) , (77) 

5 io 5 

M{rYC{r)C\r)\'^ = j^M{rfdr - ^ j\{rfC{rfdr - ^ (4/io + /i, - 3/i,) . (78) 



Using the boundary conditions for the metric functions (|3TD and (]46|) , and taking the difference 
of the eqs. (|7^ and ( |75D then establishes the first part of eq. (pB]). To prove the second 
part of (^) one starts right from the general expressions of the stress-energy components e^, 
relations (|^)-(|2|): 



00 



^^o-^^e= I drM{rfC{rf (e, - eo) = / drM{rfC{rf { + + ^^"^'^^"^ 



(79) 
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Multiplying the equation of motion for the gauge field (|TU]) by K{r) and substituting the 
J{rYK{rY term gives 



, ^4 I rw. ^2 I ' K{rf [M {rY K' {r)\ K {r) 



+ 



M(r) 



^0) 



Integration by parts in the last term of the above equation leads to 

'l-K{rf 



drM{r) 



C(r) 



+ K{r)M{r)^K'{r)\^ 



(81) 



which together with the behavior of the gauge field at the origin -^''(O) = (see (|39D) finishes 
the proof of relation (^81). 

The proof of relation (^) is simply obtained by rewriting (^) with vanishing left hand side. 

To establish ( |60|) we start directly from the definitions of the stress-energy tensor compo- 
nents, relations (B3)-(|25|): 



yUo + /Up + 2/ie 



drM{ryC{r) 



^-(>f+"<™%a(J(r)'-l)- 



12) 



Collecting derivatives in the equation of motion for the scalar field (p!8| ) and multiplying by 
J(r) gives 



J{r) - aJ{rf {j{rf - l) . 



Eliminating now the second term in the equation (^) leads to 

[M{rYC{rfj'{r)\ 



^3) 



drM{rfC{r) 



J'irf + 



J(r) + a (1 - J{rf) 



If we now expand and integrate the second term by parts we are left with 

/■oo 

/io + /ip + 2/ie = a / dr (1 - J{rf) M{rYC{rf + M{rYC{rfj{r)J\r)\'^ , (; 
Jo 

which reduces to (pO| ) when the boundary conditions for the metric ( |3TD and (|46| ) are used. 



M) 



;5) 



B Numerics 



As already pointed out, the numerical problem encountered is to find those solutions to the sys- 
tem of differential equations (|T5D-([T9D and boundary conditions for which the integral defining 
the 4-dimensional Planck-scale (|3^ ) is finite. This is a two point boundary value problem on 
the interval r = [0, oo) depending on three independent parameters (a,/3,7). Independently of 
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the numerical method employed, the system of equations ([T3|)-(]TP|) was rewritten in a different 
way in order for the integration to be as stable as possible. By introducing the derivatives of 
the unknown functions J{r), K{r), M{r) and C{r) as new dependent variables one obtains a 
system of ordinary first order equations. In the case of M(r) it has proven to be convenient to 
define = M'{r)/M{r) as a new unknown function (rather than M'(r)) since the boundary 
condition for at infinity then simply reads 

lim ^7 = - - . (86) 

r— >oo Z 

With the following definitions we give the form of the equations which is at the base of several 
numerical methods employed: 

1/5, 
2/6, 
2/3 2/7, 
2/8, 

-2 (21/7 + l) 2/5 + 2 ^ + ayi {yf - 1) , (87) 

- 42/62/7 + 2/i2/2 , 

-42/?-22/7|-f (27 + eo-e,-2e,) , 



2/1 ( 




= Ar), 


y'l 


2/2 ( 


V) 


= K{r), 


y'2 


2/3 ( 


[r) 


= M{r), 


y'3 


2/4 ( 




= ^r), 


y'. 


2/5! 


V) 


= J'{r), 


y'. 


2/6 ( 


V) 


= K'{r), 


y'. 


2/7! 


V) 


= M'{r)/M{r), 


y'l 


2/8 ( 


V) 


= AO, 


y's 



-t- 42/72/8 + ^ - ^ (27 - 4eo - e, + 3e. 



with 



.„-.,-2.„ = |(;,J-lf-(i_il£-2|, (88) 

^ 2/4 2/4 

-46o-6, + 36, = ^(2/?-l)%4^1^ + 5^ + 34. (89) 

2 ' ' 2/1 2/12/1 

This is an autonomous ordinary system of coupled differential equations depending on the 
parameters (a, /?, 7). The boundary conditions are 

yi(0) = 0, ^2(0) = 1, ^3(0) = 1, ?/4(0) = 0, 

lim|/i(r) = 1, lim|/2(r) = 0, ^^(0) = 0, y^ = 1. ^^"^^ 



In order to find solutions with the desired metric asymptotics at infinity it is useful to define 
either one (or more) of the parameters {a, f5, 7) or the constants J'{0) and K"{0) as additional 
dependent variables, e.g. yg^r) = a with y'g{r) = 0, see Before discussing the different 

methods that were used we give some common numerical problems encountered. 



• Technically it is impossible to integrate to infinity. The possibility of compactifying 
the independent variable r was not believed to simplify the numerics. Therefore the 
integration has to be stopped at some upper value of r = r,„^^. For most of the solutions 
this was about 20 ~ 30. 
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• Forward integration with arbitrary but fixed values of (a, /3, 7, J'(0) , K"{0)) turned out 
to be very unstable. This means that even before some integration routine (e.g. the 
Runge-Kutta method [^) reached r^^x, the values of some yi went out of range which 
was due to the presence of terms or quadratic and cubic (positive coefficient) terms 
in Vi- 

• Some right hand sides of (^) contain terms singular at the origin such that their sums 
remain regular. Starting the integration at r = is therefore impossible. To overcome this 
problem the solution in terms of the power series (PUD-(P^) was used within the interval 
[0, e] (for e = 0.01). 

• The solutions are extremely sensitive to initial conditions which made it unavoidable 
to pass from single precision to double precision (from about 7 to about 15 significant 
digits). However this didn't completely solve the problem. Even giving initial conditions 
(corresponding to a gravity localizing solution) at the origin to machine precision is in 
general not sufficient to obtain satisfactory precision at r^^^ in a single Runge-Kutta 
forward integration step from e to r^^^^. 



The method used for exploiting the parameter space {a, (3, 7) for gravity localizing solutions 



was a generalized version of the shooting method called the multiple shooting method |41], |42 



Hi) mi- In the shooting method a boundary value problem is solved by combining a root- 



finding method (e.g. Newton's method ^^) with forward integration. In order to start the 
integration at one boundary, the root finding routine specifies particular values for the so-called 
shooting parameters and compares the results of the integration with the boundary conditions 
at the other boundary. This method obviously fails whenever the "initial guess" for the shooting 
parameters is too far from a solution such that the forward integration does not reach the second 
boundary. For this reason the multiple shooting method was used in which the interval [e , r^^J^ 
was divided into an variable number of sub-intervals in each of which the shooting method was 
applied. This of course drastically increased the number of shooting parameters and as a result 
the amount of computing time. However, it resolved two problems: 



• Since the shooting parameters were specified in all sub-intervals the precision of the 
parameter values at the origin was no longer crucial for obtaining high precision solutions. 

• Out of range errors can be avoided by augmenting the number of shooting intervals (at 
the cost of increasing computing time). 



Despite all these advantages of multiple shooting, a first combination of parameter values 
(a, P, 7, J'(0) , K"{0)) leading to gravity localization could not be found by this method, since 
convergence depends strongly on how close initial shooting parameters are to a real solution. 
This first solution was found by backward integration combined with the simplex method 
for finding zeros of one real function of several real variables. For a discussion of the simplex 
method see e.g. [0. This solution, corresponds to the parameter values (a = 1.429965428, P = 
3.276535576, 7 = 0.025269415). 
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Once this solution was known, it was straightforward to investigate with the multiple shoot- 
ing method which subset of the (a, (3, 7)-space leads to the desired metric asymptotics. We 
used a known solution (ai, Pi, 71) to obtain starting values for the shooting parameters of a 
closeby other solution (ai + 6a, Pi + 6P, 71 + 57). This lead in general to rapid convergence 
of the Newton-method. Nevertheless, Sa, 6P, ^7 still had to be small. By this simple but 
time-consuming operation the fine-tuning surface in parameter-space, presented in section ^.2| 
and shown in Fig. ^ was found. 
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